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Abstract. In this paper we extend earlier results regarding the effects of the lower layer 
of the ocean (below the thermocline) on the baroclinic instability within the upper layer 
(above the thermocline). We confront quasigeostrophic baroclinic instability properties 
of a 2.5-layer model with those of a 3-layer model with a very thick deep layer, which 
has been shown to predict spectral instability for basic state parameters for which the 
2.5-layer model predicts nonlinear stability. We compute and compare maximum normal- 
mode perturbation growth rates, as well as rigorous upper bounds on the nonlinear growth 
of perturbations to unstable basic states, paying particular attention to the region of 
basic state parameters where the stability properties of the 2.5- and 3-layer model differ 
substantially. We found that normal-mode perturbation growth rates in the 3-layer model 
tend to maximize in this region. We also found that the size of state space available for 
eddy-amplitude growth tends to minimize in this same region. Moreover, we found that 
for a large spread of parameter values in this region the latter size reduces to only a small 
fraction of the total enstrophy of the system, thereby allowing us to make assessments of 
the significance of the instabilities. 
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1. Introduction 

Observations indicate that most of the world oceans variability is confined 
in a thin layer limited from below by the permanent thermocline. There, 
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the density is approximately uniform in the vertical but has important 
horizontal gradients. The latter imply the existence of a considerable reser- 
voir of potential energy within this layer, stored in the isopycnals tilt and 
available to feeding baroclinic instability processes ( Gill et al., 1974 1 . Haine 
and Marshall (1998) have argued that these processes are of outmost im- 
portance for the dynamics of the upper ocean. These authors pointed out 
that baroclinic instability waves can be efficient transport agents capable 
of stopping convective processes, thereby exerting a large influence in the 
thermodynamic state of the ocean. 

Because the total depth of the ocean is much larger than that of the 
upper thermocline layer, the reduced-gravity setting has been commonly 
adopted to studying the upper ocea n baroclinic instability (|Fuk amachi et al., 1995[ 
Ripa, 1995||Young and Chen, 1995||Beron-Vera and Ripa, 1997||Ripa, 1999b| 



Olascoaga and Ripa, 1999HRipa, 1999c||Ripa, 1999al|Ripa, 2000a||Ripa, 20 00b 
Ripa, 2001 ). In this setting the active fluid layer is considered as floating on 



top of a quiescent, infinitely deep layer. Olascoaga (2001) showed, however, 
that a thick — but finite — abyssal active layer can substantially alter the sta- 
bility properties of the upper ocean for certain baroclinic zonal flows, such 
as the Atlantic North Equatorial Current (ANEC) ( |Beron-Vera and Olascoaga, 2003). 
Olascoaga (2001) considered spectral (i.e. linear, normal- mode) , formal (or 
Arnold), and nonlinear (or Lyapunov) stability (Holm et al., 1985; cf. also 
Mclntyre and Shepherd, 1987) in a 3-layer quasigeostrophic (QG) model. 
Primary attention was given to the limit of a very thick bottom layer. 
The stability results were compared with those from a reduced-gravity 
2-layer (or 2.5-layer) model flOlascoaga and Ripa, 1999 ), and assessments 
were made of the influence of the deep ocean on upper ocean baroclinic 
instability. 

To make further assessments, in this paper we turn our attention to 
baroclinic instability saturation. Employing Shepherd's (1988) method, we 
establish and confront rigorous bounds on nonlinear instability saturation 
in 2.5- and 3-layer models. This method, which builds on the existence of a 
nonlinear stability theorem, has been previously used to compute saturation 
bounds in 2.5- (Olascoaga and Ripa, 19991 and 3-layer ( Paret and Vanneste, 1 996 ) 
models. In addition to considering more general model configurations than 
in these earlier works, we focus on the size of state space available for 
the growth of eddies in the region of parameter space where the mod- 
els present discrepancies in their stability properties. Also, unlike Paret 
and Vanneste (1996), who computed numerical energy-norm bounds based 
on both Arnold's first and second theorems, we derive analytical expres- 
sions for enstrophy-norm bounds based on Arnold's first theorem. Maxi- 
mum normal-mode perturbation growth rates in the 3-layer model are also 
calculated and contrasted with those in the 2.5-layer model. 
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The reminder of the paper has the following organization. Section 2 
presents the 3-layer model, from which the 2.5-layer model follows as a 
limiting case. Normal-mode perturbation growth rates are computed in 
§ 3, along with an exposition of the main results of formal and nonlinear 
stability analyses. Nonlinear saturation bounds are then derived in § 4. We 
want to remark that the number of basic state parameters that define a 3- 
layer flow is too large to be explored in full detail. To facilitate the analysis 
we reduce in some cases this number by fixing certain parameters to values 
that can be taken as "realistic," because of being found appropriate for 
a region of the ocean mainly dominated by the ANEC, which is a good 
example of a major zonal current. Section 5 presents a discussion and the 
conclusions. Appendices A and B are reserved for mathematical details 
relating to the computation of the saturation bounds in the 3- and 2.5-layer 
models, respectively. 



2. The Layer Models 



Let x denote the horizontal position with Cartesian coordinates x (east- 
ward) and y (northward), let t be the time, and let D be an infinite (or 
periodic) zonal channel domain on the (3 plane with coasts at y = ±\W. The 
unforced, inviscid evolution equations for QG motions in a 3-layer model, 
with rigid surface and flat bottom, are given by (cf. e.g. Ripa, 1992) 

d tqi = z • V% x V^i, if = 0, (2.1a) 

where ipi, being a nonlocal function of q := (qi) T and 7 := (jf) 1 " , 
uniquely determined by 
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dx dyipi = 7^, d x ipi = 



(2.1c) 



(2.1d) 



at y = ±t;W. Here, qi(x, t), ipi(x, t), and 7^ = const, denote the QG poten- 
tial vorticity, streamfunction, and Kelvin circulation along the boundaries 
of the channel, respectively, in the top (i = 1), middle (i = 2) and bottom 
(i = 3) layers. The Coriolis parameter is represented as / = /o + (3y, the 
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Nabla operator V = (d x ,d y ), and z denotes the vertical unit vector. The 
quantities 

_n Qi Hi Ho 01 Hi H , . 

i? 2 : = yi -L £ s : = 2i n := —± r 2 := — . (2.2 

/o 2 ^ 52' #2 #3 1 ^ 

where gi is the buoyancy jump at the interface of the i-th and (i + l)-th 
layers, and H := Hi + H2 with Hi the i-th layer reference thickness. 

The 2.5-layer model follows from (|2.1j) in the limit of infinitely thick 
(ra — > 0) and quiescent (^3 — ► 0) lower layer. In the latter case, (1 + 
ri)(ri/s) 1//2 i2 and R are equal to the first (equivalent barotropic) and 
second (baroclinic) deformation radius, respectively, in the limit of weak 
internal stratification (s — > 0). 

The evolution of system (|2.1|) is constrained by the conservation of 
energy, zonal momentum, and an infinite number of vorticity-related 
Casimirs, which are given by 

£:=-3<^<&>. M:=(y qi ), C := (Qfo)) (2-3) 

(modulo Kelvin circulations along the boundaries), where Cj(-) is an arbi- 
trary function and (•) := J2i Hi f D d 2 x (•). 

3. Spectral, Formal, and Nonlinear Stability 

In this paper we deal with the stability of a basic state, i.e. equilibrium 
or steady solution of (|2.1j) , of the form 

*i = -U iy , (3.1) 

which represents a baroclinic zonal flow. Here, U{ = Yhi=i 9h Ylil=i Hte,v/ fo 
= const., where Hi(y) is the thickness of the i-th layer in the basic state, 
whereas U3 is an arbitrary constant (set here to zero with no loss of general- 
ity). The following six parameters are enough to characterize the solutions 
of the 3-layer model stability problem: 

(3.2) 

The first parameter, k, is a nondimensional wavenumber of the perturba- 
tion, where k and I are the zonal and meridional wavenumbers, respectively. 
The second parameter, s, is a nondimensional measure of the stratifica- 
tion. The third parameter, b, is a planetary Charney number, namely 
the ratio of the planetary (3 effect and the topographic effect due to the 
geostrophic slope of the upper interface. Here, U s := U\ — U2 is the velocity 
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jump at the interface (i.e. the current vertical "shear"). The fourth param- 
eter, &Tj is a topographic Charney number given by the ratio of the 
topographic (3 effects due to the geostrophic slopes of the lower and upper 
interfaces. Finally, the fifth (resp., sixth), n (resp., r 2 ), parameter is the 
aspect ratio of the upper to intermediate (resp., upper-plus- intermediate 
to lower) reference layer thicknesses. The problem considered by Olascoaga 
(2001) had r\ = 1. In turn, the 2.5-layer problem treated by Olascoaga and 
Ripa (1999), which can be recovered upon making r 2 — > and ^3 — * 0, also 
had 7*1 = 1. 

Choosing a Casimir such that 5{£ + C — aftA) = for any constant a, 
the pseudo energy momentum, 

:= (A - S)(S + C- aM) = £[6^} + \{C im 8q 2 ), (3.3) 

where ip := (ipi)^, is an exact finite- amplitude invariant, quadratic in 
the perturbation Sqi(x,t) on the basic state potential vorticity Qi(y). 
Here, the symbols A and 5 stand for total and first variations of a func- 
tional, respectively, and Ci(Qi) = J dQi (a — Ui)Y(Qi) [Y is the meridional 
coordinate of an isoline of Qi] , where 

Qx = f + (b + p)yU s /R 2 , (3.4) 
Q2 = fo + [b + pn(b T -l)]yU s /R 2 , (3.5) 
Q3 = fo + (b- p 2 nr 2 6 T ) yU s /R 2 , (3.6) 

with p := (l + ri) -1 . Arnold's (1965; 1966) method for proving formal 
stability of Qi relies upon the sign-definiteness of 7i a . For evaluating the 
latter, it is useful to make the Fourier expansion e)q = ]T] fc t q(t)e lfcx sin ly, 
which implies 7i a = 5 Y^k z(*i*) T '"'aq for certain matrix H„(k, s, b, br, n, r 2 ) 
(cf. Beron-Vera and Olascoaga, 2003, §2.2.1), so that the sign of H a is 
determined from the inspection of the elements of H Q (cf. Mu et al., 1994; 
Paret and Vanneste, 1996; Ripa, 2000a). 

In figure n the regions of the (6, 6x)-space labeled "Stable" correspond 
to basic states for which there exists a such that 7i a is positive definite 
(Arnold's first theorem). The regions labeled "Possibly Stable" are locus 
of basic states for which there exists a such that 7i a is negative definite 
(Arnold's second theorem) if the channel in which the flow is contained is 
sufficiently narrow; cf. Mu (1998) and Mu and Wu (2001) for details on 
optimality issues relating to Arnold's second theorem. The results, which 
are independent of the choice of s, are presented for t\ = 0.5, a value 
estimated for the ANEC. The r.h.s. panel in this figure corresponds to the 3- 
layer model in the limit r 2 — > 0; the l.h.s. panel corresponds to the 2.5-layer 
model. Clearly, as r 2 — > the 3-layer model stable region does not reduce to 
that of the 2.5-layer model; it also requires 5ip3 — > (Olascoaga, 2001). In 
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Figure 1. Stability/instability regions in the planetary b vs. topographic &t Char- 
ney numbers space. Dark-shaded regions are the locus of positive-definite pseudo 
energy-momentum integrals. In the blank regions a pseudo energy-momentum integral 
can be found to be negative definite if the zonal channel flow is narrow enough. In the 
light-shaded regions no pseudo energy-momentum integral can be proved to be sign 
definite. 

the regions labeled "ERF no sign-definite TL a can be found. Consequently, 
the corresponding states are always unstable either through normal-mode 
perturbations or explosive resonant interaction (ERI) (Vanneste, 1995). By 
contrast, in the 2.5-layer instability problem all basic states subject to ERI 
are spectrally unstable. Finally, nonlinear stability can be proven for all 
formally stable states. Namely, the departure from these basic states can 
be bounded at all times by a multiple of the initial distance. 

For spectral stability a perturbation is assumed to be infinitesimal 
and with the structure of a normal mode, i.e. q = eqe _lfcc * + 0(e 2 ), where 
e — > 0. Nontrivial solutions for q, which satisfies H c q = 0, require condition 
det H c = to be fulfilled. This implies the eigenvalue c(k; s, b, bp, n, r%) 
to satisfy P(c) = 0, where P(-) is a cubic characteristic polynomial (cf. 
Beron-Vera and Olascoaga, 2003, appendix B). 

Figure |21 shows the 3-layer model maximum perturbation growth rate, 
max K {ftImc}, for n = 0.5, and different values of parameters r<i and s 
in the planetary b vs. topographic &t Charney numbers space. In general, 
the maximum perturbation growth rate increases with increasing T2 and 
decreasing s. As &t increases, the maximum perturbation growth rate tends 
to achieve the largest values in the region where the 2.5-layer model is 
nonlinearly stable as a consequence of Arnold's first theorem, even for 
(realistically) small values of r<i as depicted in the bottom panels of the 
figure. 

Figure |21 shows instability regions in (k, 6T)-space for b = —0.35, a value 
estimated for the ANEC, and the same values of parameters r%, r2, and s 
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Figure 2. Maximum 3-layer model perturbation growth rate in the planetary b vs. to- 
pographic 6t Charney numbers space for a fixed value of the aspect ratio ri(= 0.5) of the 
upper to intermediate reference layer thicknesses, and various values of the aspect ratio T2 
of the upper-plus-intermediate to lower reference layer thicknesses and the stratification 
parameter s. 



as in figure |^1 The area of the region of possible wavenumbers for destabi- 
lizing perturbations in the 3-layer model increases with increasing r2 and 
decreases with decreasing s. Thus the likelihood of these instabilities, which 
are not present in the 2.5-layer model, appear to be quite limited because 
they are confined only to small bands of wavenumbers for small s and r<i- 
Yet the perturbation growth rates in these bands are not negligible even 
for very small values of according to Olascoaga (2001), who explained 
these instabilities as a result of the resonant interaction between a neutral 
mode of the 2.5-layer model instability problem and a short Rossby wave 
in the bottom layer of the 3-layer model. In the next section we will see, 
however, how the existence of nonlinearly stable states contributes to arrest 
the eddy-amplitude growth, restricting the significance of these instabilities, 
at least for certain basic state parameters. 

Let us finally turn our attention to the region of (b, &x)-space where 
the two models allow for the possibility of instability. The 2.5-layer model 
acquires its maximum perturbation growth rates for b = ^p[(l — &T) r i — 
1] = — g(l + &t) an d &t < 1 + ^f 1 = 3 as s — > 0, which corresponds 
to an exact cancellation of the planetary and topographic j3 effects, i.e. 
Qi,y + Q2,y = (Olascoaga and Ripa, 1999). This result does not hold 
for the 3-layer model because of the presence of instabilities not present in 
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Figure 3. Three-layer model instability regions in the nondimensional wavenumber k 
vs. topographic Charney number 6t space, for a fixed value of the planetary Charney 
number b(= —0.35), and the same values of the aspect ratio parameters n and ri, and 
the stratification parameter s as in figure 2. 

the 2.5-layer problem. The latter instabilities are confined to very narrow 
branches in the (re, &x)-space and were also explained by Olascoaga (2001) 
as a result of the resonant interplay of a neutral mode in the 2.5-layer 
model instability problem and a short Rossby wave in the bottom layer of 
the 3-layer model. The maximum perturbation growth rates associated with 
these instabilities are larger than those of the 2.5-layer model (Olascoaga, 
2001). In the following section we will see, however, that the fraction of the 
total enstrophy of the system available for eddy-amplitude growth can be 
much smaller in the 3-layer model than in the 2.5-layer model for certain 
unstable basic states. 

4. Upper Bounds on Instability Saturation 

When a basic state is unstable, a priori upper bounds on the finite-amplitude 
growth of the perturbation to this state can be obtained using Shepherd's 
(1988) method. This method relies upon the existence of a nonlinear stabil- 
ity theorem, and the bounds are given in terms of the "distance" between 
the unstable basic state, Qf say, and the nonlinearly stable state, Qf say, 
in the infinite-dimensional phase space. 

Let 5<^(x, t) be that part of the perturbation representing the "waves" or 
"eddies," which result upon subtracting from the perturbation its zonal (i.e. 
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Figure 4- Fraction of the total potential enstrophy available for eddy-amplitude growth 
in the 2.5- (dashed lines) and 3- (solid lines) layer models as a function of the planetary 
Charney number b, with aspect ratios n = 0.5 and V2 = 0.1 (the 2.5-layer model has 
t-2 — > 0), for different values of the topographic Charney number &t- 



along-channel) average. Let § denote the space of all possible nonlinearly 
stable basic states, let ||a|| 2 := (a 2 )/Z, where Z := {(QY) 2 ), and assume 
Qi ~ QY at t = so that Z corresponds to the total enstrophy of the 
system. According to Shepherd (1988), a rigorous enstrophy-norm upper 
bound on eddy- amplitude growth, based on Arnold's first theorem, must 
have the form 

We want to mention that bounds — not treated here — on the zonal-mean 
perturbation, 5qi(y,t) := <% — 5q[, or the total perturbation can also be 
derived as <%, and Sq^ satisfy the Pythagorean relationship ||<5q|| = 
||<5q|| 2 + ||«Jq'|| 2 (cf. Ripa, 1999c; Ripa, 2000b). 

Figure 0] shows the tightest bound on instability saturation correspond- 
ing to the 3-layer model (thick curves, cf. appendix A) as a function of the 
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planetary Charney number b, for aspect ratios n = 0.5 and ri = 0.1, and 
various values of the topographic Charney number &t- The focus on the 
small value T2 = 0.1 will allow us to make comparisons with the stability 
properties of the 2.5-layer model. The latter model's bound (cf. appendix B) 
is also plotted in the figure (dashed lines) assuming n = 0.5 and the same 
values of &t- It is important to remark that the 2.5-layer bound does not 
follow from that of the 3-layer model in the limit — > 0; it also requires 
"03 = (Olascoaga, 2001). Both the 2.5- and 3-layer model bounds are 
independent of the stratification parameter s. The 2.5-layer model bound 
curves are only present in the upper three panels of the figure because the 
2.5-layer model predicts nonlinear stability as a consequence of Arnold's 
first theorem for &t > 1 + r \ 1 = 3 and any value of b (cf. also figure 
lower-right panel). 

Vertical dashed lines in each panel of figure ^indicate the values of b for 
marginal stability, denoted by b±. In the 3-layer model, 6_ = min{— p, pr\(l 
— 6t)} and 6+ = max{pri(l — b^), p^ri^bn?}, whereas for the 2.5-layer 
model, 6_ = — p and b + = pr\{l — b^). Note that in the 2.5-layer model, 
while the b^ marginal stability value remains fixed at b ~ —0.66667, the 6+ 
moves toward smaller values as 6t increases, until it collapses with 6_ at 
&t = 3 (not shown in the figure). For b > b + and b < 6_ the basic flow in 
both the 2.5- and 3-layer models is nonlinearly stable. For 6 + < b < 6_ the 
basic flow is unstable unless the zonal channel flow is narrow enough for 
Arnold's second stability theorem to be fulfilled. The latter is not always 
true in the 3-layer model case, however, since there is a possibility that a 
spectrally stable basic state could become unstable through ERI. 

Three-layer surface-confined flows are susceptible to suffer more desta- 
bilization than 2.5-layer flows. However, the state space available, (deter- 
mined by the fraction of total potential enstrophy), for eddy-amplitude 
growth in the 3-layer model tends to be smaller than the space available 
in the 2.5-layer model, at least for certain basic state parameters. This is 
evident in the upper three panels of figure 01 There is an overall tendency 
of the 3-layer model bound to decrease as &t increases. Moreover, for a 
large set of parameters this bound reduces to only a small fraction of 
the total enstrophy of the system. In these cases, the significance of the 
associated instabilities is relative. On the other hand, there are basic state 
parameters for which this fraction is not negligible. As an example not 
shown in the figure, for b = —0.35 and &t = 2.5, which are appropriate for 
a region similar to the ANEC, the fraction of total enstrophy is about 45%, 
which is not negligible. Of course, when the upper bounds are not small 
enough, no unambiguous conclusion can be drawn about the significance of 
an instability. 

Figure ^ also shows that the 3-layer model bound can be significantly 
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small for certain potentially ERI unstable flows (cf. lower two panels in 
the figure). This also allows us to make an unambiguous assessment of 
the significance of these type of instabilities in the sense that they can be 
certainly negligible for some basic state parameters. 

Before closing this section, two points deserve additional discussion. 
First, the result that the bounds for the 3-layer model with a very thick 
deep layer are smaller than the 2.5-layer model bounds in the region of 
parameters where the two models share similar instability properties might 
seem at odds with the fact that the 3-layer model is less constrained than 
the 2.5-layer model, which allows for the development of more unstable 
states. However, we believe that this result should not be surprising inas- 
much as the space over which the minimization is carried out is larger in 
the 3-layer model than in the 2.5-layer model, which offers the possibility of 
finding tighter bounds (cf. Olascoaga, 2001). Second, Paret and Vanneste 
(1996) were not able to draw a conclusion on the significance of ERI in- 
stability as in the present paper. These authors computed energy-norm 
saturation bounds, according to both Arnold's first and second theorems, 
using numerical minimization algorithms. These bounds, whose analytical 
computation appears to be too difficult, were not found to minimize at 
basic state parameters for which ERI instability is possible. The analytical 
minimization involved in the derivation of the enstrophy-norm of this paper 
has shown that the tightest bounds are obtained using stable basic states 
whose parameters have quite spread numerical values. The minimization 
thus requires to search for a solution in a considerably large space, making 
numerical computations extremely expensive. This might explain the diffi- 
culty of Paret and Vanneste (1996) to find tighter bounds for potentially 
ERI unstable basic flows. 

5. Concluding Remarks 

A previous study showed that the quasigeostrophic baroclinic instability 
properties of a surface-confined zonal current may differ substantially be- 
tween a 2.5-layer model and a 3-layer, if the former is considered to be a 
simplified 3-layer model with a very thick deep layer. For certain basic state 
parameters, the 2.5-layer model predicts nonlinear stability whereas the 3- 
layer model spectral instability. That study thus suggested that the effects 
of the deep ocean on the baroclinic instability of the upper thermocline 
layer of the ocean may be important for certain currents. 

In this paper we have made further assessments of the importance of 
the deep ocean on upper baroclinic instability. We have achieved this by 
analyzing (i) maximum normal-mode perturbation growth rates and (ii) 
rigorous enstrophy-norm upper bounds on the growth of perturbations to 
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unstable basic states, in both 2.5- and 3-layer models of baroclinic instabil- 
ity. The new results show that instabilities, which the 3-layer model predicts 
in the region of basic state parameters where the 2.5-layer model predicts 
nonlinear stability, appear to maximize their growth rates. At the same 
time, however, the saturation bounds tend to minimize in this same region 
of basic state parameters, thereby reducing the size of state space available 
for eddy-amplitude growth. Moreover, for a large subset of parameters in 
the region, the latter reduces to only a small fraction of the total enstrophy 
of the system. In these cases we have been able to make unambiguous 
assessments of the significance of the associated instabilities in the sense 
that they can be certainly negligible. 

We close remarking that the important issue of making assessments of 
the accuracy of the saturation bounds as predictors of equilibrated eddy 
amplitudes is still largely open. This cannot be addressed without per- 
forming direct numerical simulations. The importance of this subject relies 
upon the potential of the bounds in the architecture of transient-eddy 
parametrization schemes. The treatment of these issues are reserved for 
future investigation. 
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A. Three-Layer Model Bounds 



Upon minimizing the r.h.s. of (|4.1|) over all stable states, we have been 
able to find, in addition to the trivial bound max various sets 

of possible bounds. A first set involves 9 possibilities, for which — 
max{Qf y } and is given by 



■v 



max||<5q'|| 2 = j Q lv { - Q lv +Q ^ ( A .l) 

Ylj(Qj,y) 2 ' for i ^ / = 1, 2, 3. A second set involves other 9 possibilities, 
for which max{Qf y } = Q\^ y = Q%, y an d is given by 

max||5q'f =\ ±(QF 1)2/ - Q^) 2 (A.2) 
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4 EiiQly) 2 , fOT * + h,h, where {h,I 2 } = {1, 2}, {2, 3}, {1, 3}. Another 



possibility finally results for Qf y = Q| = Qg , and is given 



max 



INI 



i - 



Ql,yQ2,y + Ql,yQs,y + Q^yQ^y 



(A.3) 



The tightest bound follows as the least continuous bound of the above 
20 possible bounds in the 4-dimensional space of unstable basic state pa- 
rameters, with coordinates (6, br, n, ^2) (the bounds are independent of 
s). 



B. Two-and-a-Half-Layer Model Bounds 

In the 2.5-layer model (r 2 — > and ^3 - ► 0) the least bound in the 3- 
dimensional space of unstable basic state parameters, with coordinates 
(b,br,ri), is given by 

-±Ql,y(Ql,y + Q2,y) ^ ~ P < b < h 
h(Qty-Qly) 2 if 61 < 6 < 6 2 (B.l) 

-4Qly(Ql,y + Qly) if h 2 < b < -r lP (b T ~ 1) 

- [(Qi, y ) 2 + (Qly) 2 ], where ^ := - \p [n(6 T - 1) + 3] and 6 2 := 

x [3ri(6r — 1) + !]• This result extends to arbitrary n that of Olascoaga 

and Ripa (1999). 



II r /l|2 

max oq = < 
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